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Abstract 

A new exact nonlinear Newtonian solution for a plane matter flow superim- 
posed on the isotropic Hubble expansion is reported. The dynamical effect of cosmic 
vacuum is taken into account. The solution describes the evolution of nonlinear 
perturbations via gravitational instability of matter and the termination of the per- 
turbation growth by anti-gravity of vacuum at the epoch of transition from matter 
domination to vacuum domination. On this basis, an 'approximate' 3D solution is 
suggested as an analog of the Zeldovich ansatz. 



1 Introduction 



Three decades ago, Zeldovich (1970) published in this Journal his now famous nonlinear 
theory of gravitational instability formulated in terms of Newtonian mechanics. The the- 
ory was first used for pancake cosmogony by Zeldovich and his group (see for a review 
Shandarinfe Zeldovich 1989). Later on it was realized that the theory can be judiciously 
applied to a wide range of cosmological scenaria; in particular, the dynamics given by 
the theory underlies gravitational clustering (Peebles 1993). The Zeldovich theory pro- 
vides also an effective analytical tool for optimization of cosmological simulations (Melott 
1993). 

The Zeldovich theory is based on an exact analytical solution which describes a plane 
pressure-free matter flow superimposed on the regular Hubble expansion: 

xi = a{t)xi + b{t)l3{xi), X2 = a{t)x2, X3 = a(t)x3, (1) 

where x = {xi,X2,X3) are the Eulerian coordinates (Cartesian) of an element ('particle') 
with a Lagrangian coordinates x = (Xi? X2, Xs)) ^ind /3(xi) is an arbitrary function of only 
one coordinate Xi- The first term in the first equation represents — together with the two 
other equations — the unperturbed isotropic solution for the Hubble flow, and the second 
term in the first equation represents a perturbation flow which depends on Xi and time t. 
The solution was obtained by Zeldovich for the parabolic expansion (a spatially flat 

/3 \2/3 

model), in which the scale factor a(t) = aoi-Hotj , where and Hq are the present 

values of a(t) and Hubble constant, b{t) =t^/^. The perturbation corresponds to the grow- 
ing mode of Lifshits (1946) linear solution for gravitational instability of non-relativistic 
matter and has the same time dependence as in the linear theory. 

A decreasing mode of Lifshits theory can also be introduced to the solution (Zentsova 
and Chernin 1980), and then one has instead of the first equation in Eq.(|I]): 

xi(xi, t) = a{t)xi + b{t)f3ixi) + c(t)7(Xi), (2) 

where c{t) = t~^^^, and 7(xi) is the second arbitrary function of Xi- For this more general 
solution, which is also nonlinear and exact, the density is 



dxi dxi 



(3) 



In Zeldovich (1970) work, a 3D 'approximation' is suggested, in which 3D perturbation 
flow is assumed to have the same structure and time dependence as the plane perturbation 
flow in the exact solution of Eq.(|l]). This ansatz has been justified by many numerical 
simulations (see again Shandarin and Zeldovich 1989). The corresponding exact analytical 
3D nonlinear solution of the problem has not been found yet and there is little hope to 
find it. 

In this Letter, we report a generalization of Zeldovich theory by accounting for the 
dynamical effect of the cosmological constant A, or cosmic vacuum. We give a new exact 
nonlinear Newtonian solution with the same symmetry as in Eqs.(|l|-^, but for a non- 

Ac2 

zero cosmological constant, or vacuum density p\ = . The observational basis for 
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this approach is due to the recent studies of distant type la supernovae (Perlmutter et 
al. 1998, Riess et al. 1999) which indicate that the vacuum density (in the units of the 
critical density) is Qa = 0.7 ±0.1, which is in concordance with the bulk of observational 
data that come also from the cosmic age, the cosmic microwave background anisotropy 
in combination with cluster dynamics, etc. (Carrol et al. 2000, Wang et al. 2000). 

In spherical symmetry, the nonlinear solution of the same problem (without vacuum) 
is given in the frame of the well-known Tolman-Bondi-Lemaitre model (see Lahav et al. 
1991 and the references therein). 



2 Basic Equations 

Following Zeldovich (1970), we consider a plane perturbation flow of pressure-free matter 
and search for a solution in the Lagrangian form 

xi = a{t)xi + S{xi, t), X2 = a{t)x2, = a(t)x?,. (4) 

The equation of motion for the 'unperturbed' function a{t) in the background Fried- 
man solution contains the matter and vacuum densities 

a = —pca, (5) 



where the effective gravitating density 

Pg = Po^ +Pa + 3— = po^ - 2pA- (6) 



Here po is the present-day mean matter density. Pa is the vacuum pressure, and the 
equation of state of vacuum is Pa = ^Pac^- 
The first integral of Eq.(^, 



a'=^-^L§+PA]a'-kc\ (7) 



is the Friedman cosmology equation, where /c = 1, 0, —1 for elliptic, parabolic and hyper- 
bolic dynamics, respectively. 

The perturbation plane flow is considered in the Newtonian approximation. The basic 
equations for the flow of Eq.(4) are the Euler equation of motion, the continuity equation 
and the Poisson equation: 

V + {vV)v = -V<^, (8) 
P + V(p«)=0, (9) 
= 4nG{p-2pA). (10) 

Here v and p are perturbed velocity and density. The dynamical effect of vacuum on 
the flow is taken into account in the Poisson equation by the vacuum density pA which 
enters also the Eq.(KD. The velocity component Vi depends on Xi only. 



Now we use the Lagrangian coordinates instead of the Eulerian coordinates in accor- 
dance with the relations (prime denotes derivative on Xi) 
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X a + S' dxi a 
As a result the continuity equation can be presented in the form 

p a + 6' 



(11) 
(12) 

(13) 



+ 

p a + o' 



0. 



Its solution is 



P 



a'^ a + S' 



Similarly Eq.(^ is reduced to the form 

Substituting this relation into the Poisson equation gives 

-^-2- = 47rG(p-2pA). 
a + a 



(14) 



(15) 



(16) 



(17) 



With the use of Eqs.(^,(|^ and (|l5l), we obtain from Eq.([T7|) rather simple relation 
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S' 



(18) 



Assuming, as usually, that this relation holds for both S' and S, we cross then from the 
derivatives with respect to t to the derivatives with respect to s = a{t)/ao and finally 
obtain the equation that controls the time behaviour of the flow: 



Po + Pas'" 



kc^ 3 



a'o SnG J ds^ 



(19) 



Eq.([T9|) can be solved numerically with the use of a standard procedure. It is most 
important that for the parabolic {k = 0) motion, an exact solution of the equation can be 
found in an explicit analytical form. 



3 Exact solution 



Let us denote q = Pa/ Po = ^a/^o where and fio are the present vacuum and matter 
densities in the units of critical density. For k = one has J^a + = 1 and 



l + q 



1 



(20) 



Introducing u = qs^ = q[a{t) / Qq]'^ , where now 

a{t) = aoq'^'^u^'^it), u{t) = sinh^ 



-at 



(21) 



and a = y/Tl^Ho = \ q/{l + q)HQ = JSttGpa/S, we transform Eq.([l9D to the form 



d6_ 

du 



{l+u)5 = 0, 



(22) 



Eq.(p2|) can be reduced to hypergeometric equation, and its general solution is given 
in terms of hypergeometric functions. Accordingly, one has the following solution: 



6 = B{qs')h{xi)+C{qs')fc{xi] 



(23) 



The solution contains two modes of the perturbation flow (we call them 1 and 2) and 
includes /b and fc as arbitrary dimensionless functions of the Lagrangian coordinate Xi- 
The second mode is given by an elementary function: 



C{u) 



VT+u 



u 



1/6 



(24) 



As for the first mode, B{u), it is given in terms of hypergeometric functions: 
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u 



2' 6 1 + M 



CoC{u) - -D{u) 



if < M < 1, 
if M ~ 1, 
if u > 1, 



(25) 



where Cq = —=T — T 77 and 
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D(u) 



u 



1/3- 



u 



- [coC(m) - B{u) 
5 



(26) 



Fig.l demonstrates the behaviour of i?, C and D as func- 
tions of time variable u (double logarithmic scale) in the most interesting "transition" 
range of the argument u, e.g. in the era, when the evolution of the flow transits from the 
initial epoch of dynamical domination of matter gravity to the later epoch of dynamical 
domination of vacuum anti-gravity. 



IgB, IgC, IgD 
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Fig.l. The functions B{u), C{u) and D{u) 

It is also easy to see from Fig.l the dependence of the functions on the redshift z, since 
u= g(a/ao)^ = g(l + z)~^. 

Simple asymptotic formulae for the three functions at the epochs before and after the 
transition era are given in Table 1. In this table the asymptotic dependence of a and u 
on time is also shown. 

Table 1. Asymptotic behaviour of a(t), u{t), B{u), C{u), D{u) 





t 


t ^ oo 


a{t) 

u{t) 

B{u) 

C{u) 

D{u) 


aog-i/3(3at/2)2/3 
(3at/2)2 

^2/3 
^-1/6 

4coM-^/V5 


aog-^/^2-2/3e"* 
2-2exp(3at) 

„l/3 
^-1/3 



4 Discussion 

1. The exact solution ( p^) contains the Zeldovich solution of Eq.(l) and the exact solution 
of Eqs.(2,3) as particular cases for pA = 0. It is also seen from Table 1 that Eqs.(l,2) 
represent an asymptotic of our solution in the limit t — ^ 0, a{t) 0, when the dynamical 
effect of vacuum can be neglected. Indeed, in this limit, the two modes of our solution 
vary as 

B{uit)) oc b{t) = t^'\ C{u{t)) oc c(t) = t-^l\ (27) 

It may be easily found with this asymptotic that the arbitrary functions of the solution 
of Eqs.(l,2) are related to the arbitrary functions of our solution in a simple way: 



/3 \^/^ /3 V^/^ 

/?(Xi)=(^2"j /b(Xi)> 7(Xi) = (^2"j /c(Xi)- 



(28) 



2. In the opposite asymptotic region, in the limit when t — oo, a(t) — oo, vacuum 
dominates; since the gravity of matter can be then neglected, matter may be considered as 
a gas of test particles moving freely on the vacuum background. In this limiting case, the 
time behaviour of the perturbation flow differs drastically from the Zeldovich solution. 
Indeed, for large times one has from our solution: 



S oc B{u) (xC{u) (X. a{t) oc e' 



at 



(29) 



As we see, the time behaviour of the perturbation flow is exactly the same as that of the 
unperturbed Hubble flow, in this limit. Accordingly, the whole nonlinear (and arbitrary 
non-uniform in density — see below) flow expands in a regular isotropic manner with the 
Hubble linear law. 

3. Our solution allows a special case, in which functions /b and /c are related to each 
other as /b = — /c/co(= — 4/d/5); the asymptotic of this special solution at a{t) — >^ oo is 
different from both Eqs.(l,2) and the equation above: 



5 oc D{u) oc a{t) ^ 



oc e 



-at 



(30) 



The special solution describes the asymptotic "adiabatic cooling" of the nonlinear flow. 

4. According to Eq.(^), the perturbation of the matter density relative to the per- 

3 p (5' 

turbed density e = — = . One may consider also the density perturbation relative to 

p a 

the unperturbed density: 



Po 



6' 



a + S' a + S' 



e. 



(31) 



The asymptotics of the relative density perturbation e^,, when t ^ and t — oo, are 
given in Table 2 for the two modes of the general solution (denoted in the Table as 1,2) 
and for the special case 3. 



Table 2. Asymptotics of 



Mode 


t - 


> 




t ^ OO 


1 


1/3 




2/3 




-/b 

ao 


^^ao + /^cogV3 


2 




1 




1/3 

/c 


3 




1 




1/3 

^2^/3 ^-2at 

ao 



If one considers e at small times, it may be seen that its behaviour is the same as 

g^/^ /3 

for the first mode of the solution. For the the second mode e ~ — I -at ] . For 



the third mode we must change fc to /d and add factor 4co/5. 



ao 



When t oo then e coincides with for mode 3 and does not contain the second 
terms in denominators of modes 1,2. 

5. In the velocity perturbation, there is only one component, 6vi = 6, and the relative 
velocity perturbation is 

Svi 1 86 
Xia Xi9a 

The asymptotics of relative perturbations of velocity for the same three cases 1,2,3, 
as in Table 2, are given in Table 3. 



(32) 



Table 3. Asymptotics of 5*f i 



Mode 


t 


t — > oo 


1 




Xl 


2 




Xl 


3 


15 Xl 


_ /d 24/3^1/3g~2Qt 
Xl 



(We omit the argument Xi of functions / in Table 3 and of their derivatives in Table 2.) 

This analysis shows that gravitational instability is terminated and the nonlinear 
perturbations, that are growing at earlier epoch, are then freezing out or even being 
"adiabatically cooled", when vacuum starts to dominate dynamically; this effect is clearly 
due to anti-gravity of vacuum. 

6. It is worth to find a relation between the nonlinear solution and the behaviour 
of linear perturbations on the vacuum background. The general solution for linear per- 
turbations may easily be found with the method suggested by Zeldovich (1965) for the 
Lifshits-type perturbations (see for comparison Heath 1977, Lahav et al. 1991): 



8 = F^{Xi)a{t) 



dt 



+ F2(xi)a(t), 



(33) 



where Fi, F2 are arbitrary functions of Xi- Two these linearly independent linear solutions 
are directly related to modes 1,2 of the nonlinear solution obtained here if we use a{t) 
from (121]): 



a{t) 

For Mode 3 one has: 



dt 2 q 



1/3 



d^(t) 5 aoo;^ 



B{u{t)),a{t) = aoq'^/^aC{u{t)). 



D{u{t)) = 2c?aQq'^l^a 



dt 



(34) 



(35) 



7. Finally, on the basis of our solution, a 3D 'approximate' (in Zeldovich's (1970) 
sense) solution may be suggested as a generalization of the Zeldovich ansatz by taking 
into account the dynamical effect of cosmic vacuum: 



(36) 



where the time functions B and C are given by our nonlinear solution, are arbitrary 

vector functions of the Lagrangian 3D coordinates X- 

We expect that this new generalized 3D approximation can effectively be used for 
the theories of the large-scale structure formation in ACDM cosmologies, - in the same 
way as the Zeldovich original ansatz has been used for cosmologies with zero cosmological 
constant. 

It will be interesting to investigate the influence of radiation on the evolution of 
perturbations. 
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